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3  INTRODUCTION 

This  report  summarizes  the  work  performed  at  von  Karman  Institute  to  test  the 
UHTC  samples  in  the  VKI  Plasmatron  Facility.  A  DVD  attached  to  this  report  completes 
the  information  given  in  this  report. 


4  EXPERIMENTAL  SETUP 

4.1  Facility 


Figure  1:  The  VKI  Plasmatron  Facility 


The  complete  VKI  Plasmatron  Facility  [1],  [2],  [3]  is  sketched  in  Fig.  1.  It  is 
configured  with  an  ICP  torch,  of  160  mm  diameter  suited  for  material  testing.  The  torch  is 
mounted  inside  a  support  enclosure,  which  is  fixed  on  a  side  of  the  test  chamber.  The  test 
chamber  is  a  2.5  m  long,  1.4  m  diameter  vessel  equipped  with  multiple  portholes  and 
windows  to  allow  maximal  flexibility  and  unrestrained  optical  access  for  plasma 
diagnostic  techniques.  Inside  the  enclosure,  the  samples  and  probes  are  mounted  on  a  fast- 
injection  system. 
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The  plasma  jet  is  collected  at  the  outlet  of  the  enclosure  and  cooled  in  the  heat 
exchanger  to  a  maximum  temperature  of  50°C  to  protect  the  vacuum  plant  from 
overheating  damage.  The  vacuum  plant  consists  of  three  volumetric  vacuum  pumps, 
which  allow  operating  pressures  between  1  hPa  and  atmospheric  pressure  with  a 
maximum  flow  rate  of  3000  m3/h.  A  Roots  pump  can  be  inserted  in  the  circuit  to  achieve 
lower  operating  pressures.  Exhaust  gases  are  then  vented  to  the  atmosphere  through  a 
stack. 

The  Plasmatron  is  equipped  with  a  1 .2  MW,  400  kHz,  high-frequency  generator  of 
relatively  new  solid-state  technology,  using  thyristors  and  MOS  inverters  instead  of 
vacuum  tubes.  An  extensive  closed  circuit  cooling  system  using  de-ionised  water  protects 
all  facility  parts  from  melting  due  to  the  plasma  heat,  which  is  dissipated  through  three  dry 
air  coolers  located  on  the  roof  of  the  installation.  For  air  plasma  operation,  the  facility  is 
connected  to  the  VKI  compressed  air  supply. 

4.2  Instrumentation 
4.2.1  Heat  Flux  Probe 

This  probe  (Fig.  2)  is  made  mainly  of  copper  and  its  external  shape  is  the  same  as 
the  sample  holder  to  conserve  all  the  parameters  of  similitude  (dynamic  pressure, 
stagnation  heat  flux  and  velocity  gradient).  The  walls  are  cooled  by  water  and  a 
calorimeter  described  in  the  next  part  is  placed  in  the  centre  of  the  front  face.  The  heat 
flux  values  are  detennined  from  the  water-cooled  calorimeter  as  described  in  section 
4.2.1. 1. 

In  a  plasma  heating  experiment,  there  are  three  possible  contributions  to  heat  flux 
for  a  test  sample  placed  in  the  stagnation  point  configuration.  These  are  convective, 
chemical  (surface-catalyzed  recombination)  and  radiative  heating.  Note  that  the 
calorimeter  has  a  polished  copper  surface,  which  is  intended  to  improve  the  sensitivity  to 
chemical  heat  flux  rather  than  radiative  heat  flux. 

Note  that  the  heat  flux  is  measured  on  a  cold  wall  made  of  copper,  which  is 
considered  to  be  fully  catalytic  for  any  plasma  (air  or  C02),  and  it  is  polished  before  each 
application  to  give  a  reference  status  for  the  surface. 


Figure  2:  Sketch  of  the  heat  flux  probe,  with  water-cooled  calorimeter  installed 

inside 
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4.2. 1.1  Copper  Water-Cooled  Calorimeter 

The  calorimeter  is  used  to  measure  the  heat  flux  at  the  stagnation  point  of  the  heat 
flux  probe  (Figure  3).  It  measures  the  fully  catalytic  heat  flux  on  a  cold  wall  («300K)  by 
measuring  the  mass  flow  and  the  temperature  difference  in  the  cooling  water  supply  and 
return  lines. 

The  water  temperature  is  measured  in  the  inlet  and  the  outlet  of  the  calorimeter  by 
using  two  type  K  thennocouples  and  the  mass  flow  of  the  cooling  water  is  controlled  by  a 
rotameter  and  is  measured  before  every  test.  In  order  to  have  correct  measurement,  the 
sidewalls  of  the  calorimeter  have  to  be  adiabatic.  An  insulator  is  installed  between  the 
calorimeter  and  the  wall  of  the  probe  for  this  purpose.  Then  the  heat  flux  is  given  by: 

q  s 

where 

in  :  mass  flow  in  the  calorimeter  [kg/s] 

Cp  :  heat  capacity  of  water  [J/kg.K] 

Tout :  temperature  at  the  outlet  of  the  calorimeter  [K] 

Tin  :  temperature  at  the  inlet  of  the  calorimeter  [K] 

S:  front  surface  area  of  the  calorimeter  [m  ] 

A  large  number  of  tests  have  been  carried  out  with  this  calorimeter  at  VKI,  and  it 
has  proven  to  be  both  robust  and  reliable. 


K 


Water  in 

n 

i 

a 

1 

/  Water  out 


Figure  3:  Water-cooled  calorimeter 


4.2.2  Sample  Holder 

The  Sample  holder  (Fig.  4)  is  composed  of  two  coaxial  tubes  in  which  water 
circulates  for  cooling.  At  the  extremity,  a  support  for  the  sample  made  of  SiC  (Fig.  5)  is 
attached  to  the  holder  body  by  three  metallic  pins.  The  sample  is  installed  in  this  support 
and  held  with  a  second  set  of  three  pins.  An  insulator  (kapirok)  is  placed  between  the 
sample  and  the  sample  holder. 
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4.2.3  Pyrometer 

A  two-color  pyrometer  (Raytek  Marathon  Series  MR1SC,  Serial  #:  46667,  Lot  #: 
2360430101,  Temperature  Range:  1000-3000°C)  is  used  to  measure  the  surface 
temperature  of  the  sample  during  the  test.  This  device  measures  the  IR  emission  from  the 
surface  at  two  different  wavelengths  that  are  fairly  close  to  each  other.  If  the  emissivity,  s, 
of  the  material  is  considered  constant  over  the  two  wavelengths  an  absolute  temperature 
measurement  can  be  determined  from  a  Planck’s  law  analysis  for  a  grey  gas.  The 
acquisition  speed  is  1  Hz.  The  wavelength  range  of  the  pyrometers  is  0.8  to  1.1  microns 
for  the  first  detector  (used  for  one  colour  measurement)  and  1.0  to  1.1  microns  for  the 
second  detector. 

The  pyrometer  is  calibrated  using  a  black-body  radiation  source  to  take  into 
account  the  effect  of  the  window.  The  resulting  calibration-curve  is  as  follows,  with  an 
uncertainty  of  +  10C: 

Treal[C]  =  0.887Tpyrometer[C]  +  65.687 

The  angle  between  the  sample  and  the  pyrometer  is  20  degrees  in  horizontal  plane 
(Angle  ‘b’  in  Fig.  6)  and  22  degrees  in  vertical  plane  (Angle  ‘a’  in  Fig.  6). 
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4.2.4  Infrared  Camera 

A  FLIR  infrared  camera  is  used  at  the  range  of  500C  ~  1300C  to  record  the  surface 
temperature  of  the  samples.  A  Zn-Se  window  installed  at  Plasmatron  wall  allows  the 
infrared  radiaton  to  pass.  The  infrared  camera  is  not  calibrated  and  the  results  are 
emissivity  dependant.  The  wavelength  range  of  the  infrared  camera  is  8  to  9  microns. 

The  angle  between  the  sample  and  the  infrared  camera  is  26  degrees  in  horizontal 
plane  (Angle  ‘b’  in  Fig.  6)  and  33  degrees  in  vertical  plane  (Angle  ‘a’  in  Fig.  6). 


4.2.5  LIF  Diagnostics 

Early  in  the  development  of  plasma  wind  tunnels  for  aerospace  material  testing 
and  development,  attempts  were  made  to  diagnose  the  plasma  state  using  emission 
spectroscopy.  Applications  of  emission  spectroscopy  continue  in  plasma  facilities, 
however,  its  implementation  and  the  subsequent  interpretation  of  the  spectral  intensities  in 
terms  of  thermodynamic  and  chemical  parameters  of  the  plasma  has  not  been  universally 
successful.  A  number  of  applications  have  demonstrated  the  need  for  acquiring  radial 
profile  measurements,  owing  to  the  influence  of  thermal  and  chemical  gradients  along  the 
line  of  sight.  To  overcome  this  limitation,  a  number  of  research  groups  have  applied 
spatially  resolved  spectroscopic  techniques  to  characterize  plasma  flows  over  a  wide 
range  of  flow  speeds  and  test  gases.  Laser-induced-fluorescence  (LIF)  is  one  such 
approach  that  has  shown  promising  results  in  a  number  of  plasma  flow  characterization 
applications.  LIF  produces  spatially  resolved  measurement  at  the  location  of  the 
intersection  of  the  laser  beam  and  the  fluorescence  detection  optical  system.  Furthermore, 
since  fluorescence  is  used  to  monitor  the  excitation  process,  the  information  about  a 
particular  flow  species  obtained  from  the  LIF  signals  is  representative  of  the  ground 
electronic  state  population  rather  than  an  excited  state  population,  which  may  not  be  in 
Boltzmann  equilibrium  with  the  ground  state.  For  applications  involving  air  plasmas, 
typical  target  species  include  NO,  N  and  O  [9]. 
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5  TEST  CONDITIONS  AND  FACILITY  CALIBRATION 

The  purpose  of  the  test  campaign  is  to  keep  the  samples  under  plasma  flow  at  steady 
state  conditions  based  on  facility  power  and  sample  front  face  temperature.  Some  initial 
tests  are  performed  using  the  same  sample  to  see  which  power  corresponds  to  which 
surface  temperature  and  also  to  detect  the  problems,  if  any. 

Three  different  samples  are  used  for  this  purpose,  which  are  samples  no  17,  18  and 
19.  All  tests  are  performed  at  50mbar  static  pressure  conditions.  Sample  17  is  exposed  to 
air  plasma  flow  at  a  mass  flow  rate  of  16  g/s  and  at  a  power  level  varying  between  150kW 
and  200kW.  The  surface  temperature  values  have  varied  between  1260C  and  1400C.  The 
same  sample  is  afterwards  exposed  to  a  plasma  flow  of  8g/s  mass  flow  rate  at  similar 
power  conditions,  however  the  SiC  cover  is  broken  during  these  tests.  The  trials  with  the 
other  samples  (no  18  and  19)  resulted  similarly;  the  SiC  cover  is  broken  eventually 
although  a  steady  state  front  face  temperature  of  1460C  is  observed  at  175kW  power. 
Sample  19  is  also  tested  with  a  longer  SiC  cover  and  at  higher  pressure  (lOOmbar)  and 
higher  mass  flow  rate  (16g/s).  The  SiC  cover  is  not  broken  during  this  test. 

The  results  of  the  calibration  phase  forced  us  to  change  the  conditions.  For  the 
actual  tests,  the  following  conditions  are  respected: 

■  A  longer  SiC  cover  is  used 

■  The  static  pressure  is  set  to  lOOmbar 

■  Mass  flow  rate  of  air  is  set  to  16g/s 

■  Samples  are  injected  at  low  power  (~1  lOkW)  and  the  power  is  increased  slowly  to 
its  target  value  in  ~1  minute. 

■  The  power  is  decreased  slowly  at  the  end  of  the  test  to  ~100kW,  till  the  surface 
temperature  decreases  below  1000C. 


6  EXPERIMENTAL  RESULTS 

Table  1  presents  a  summary  of  all  the  tests  perfonned.  As  a  precaution,  the  first  tests 
are  started  at  lower  power  (150kW)  and  the  power  is  increased  as  the  tests  are  successful. 
Some  conditions  are  repeated  by  changing  the  total  duration,  where  the  purpose  is  to  see 
the  effect  of  duration.  No  significant  mass  change  has  been  observed  for  any  of  the 
samples  tested. 

The  heat  flux  is  determined  by  a  water-cooled  calorimeter  and  the  surface 
temperature  is  detennined  by  the  two-colour  pyrometer.  The  temperature  values  given  in 
Table  1  are  already  calibrated  values.  All  the  tests  are  carried  out  without  any  problems 
except  for  the  last  test  (Test  no  16,  Sample  4),  where  the  Plasmatron  facility  has  stopped 
after  15  minutes  of  testing.  Since  the  purpose  of  this  test  was  to  keep  the  sample  in  plasma 
flow  for  20  minutes,  the  facility  is  restarted  and  the  sample  is  kept  under  plasma  flow  for 
another  5  minutes. 
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Table  1  Summary  of  Tested  Samples 


Test 

no 

Sample 

Date 

Time 

pstat 

[mbar] 

Power 

[kW] 

H.  Flux 
[kW/m2] 

p  dyn 
[Pa] 

Duration 

Tsteady 

[C] 

5 

20 

31/5/2006 

14:30 

100 

150 

510 

24 

10  mins 

1241 

6 

16 

31/5/2006 

16:00 

100 

160 

715 

25 

10  mins 

1387 

7 

15 

1/6/2006 

10:05 

100 

160 

730 

26 

5  mins 

1356 

8 

14 

1/6/2006 

11:25 

100 

160 

755 

26 

20  mins 

1387 

9 

13 

1/6/2006 

14:00 

100 

170 

905 

28 

10  mins 

1458 

10 

12 

1/6/2006 

15:35 

100 

170 

875 

27 

20  mins 

1436 

11 

11 

2/6/2006 

11:00 

100 

170 

900 

27 

5  mins 

1456 

12 

5 

6/6/2006 

15:45 

100 

160 

710 

24 

16  mins 

1352 

13 

6 

7/6/2006 

11:15 

100 

190 

1095 

31 

20  mins 

1512 

14 

10 

7/6/2006 

14:25 

100 

190 

1050 

31 

10  mins 

1503 

15 

9 

7/6/2006 

16:05 

100 

210 

1225 

38 

10  mins 

1574 

16 

4 

8/6/2006 

12:05 

100 

210 

1165 

39 

15+5  mins 

1547 

The  temperature  evolution  of  samples  12  and  6  (tests  10  and  13)  are  presented  as 
examples  in  Figures  7  and  8,  respectively.  As  it  can  be  seen,  the  maximum  fluctuations  in 
temperature  are  in  the  order  of  +15K.  As  an  example  of  repeatability,  the  temperature 
evolution  of  sample  4  (Test  16)  is  presented  in  Figure  9.  The  Plasmatron  had  stopped  after 
15  minutes  of  testing  of  this  sample,  and  then  it  was  restarted  for  an  additional  5  minutes. 
In  Figure  9,  it  can  be  observed  that  the  steady  state  temperature  reaches  to  its  original 
value  after  the  sample  cools  down  and  the  Plasmatron  is  restarted. 

Similar  plots  for  all  the  samples  tested  can  be  found  in  the  DVD  attached  to  this 

report. 


Sample  12  -20  mins 


time  [sec] 


Figure  7.  Temperature  Evolution,  Test  10,  170kW power 
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Figure  8.  Temperature  Evolution,  Test  13,  190kW power 


Sample  4  - 15+5  mins 

1900 
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Figure  9.  Temperature  Evolution,  Test  16,  210kW power 

The  infrared  camera  is  employed  during  the  tests  13,  14  and  15  (Samples  6,  10  and 
9)  in  addition  to  the  pyrometer.  The  purpose  of  using  the  infrared  camera  is  to  detect  if  the 
emissivity  is  changing  in  time.  For  this  purpose,  measurements  are  taken  every  two 
minutes  and  a  “representative  emissivity  value”  is  obtained  from  the  software  of  the 
infrared  camera  such  that  the  infrared  camera  temperature  value  is  equal  to  the  two-colour 
pyrometer  temperature  value.  It  should  be  noted  that  this  emissivity  value  should  not  be 
taken  as  the  absolute  emissivity  value  of  the  material.  The  results  are  presented  in  Table  2. 
A  slight  increase  in  emissivity  can  be  observed  at  the  beginning  of  the  test,  but  then  the 
emissivities  of  the  materials  remain  constant. 

During  testing,  green  emission  is  observed  from  samples  once  surface  temperature 
has  reached  a  high  value.  Emission  gradually  disappears,  and  it  is  thought  to  be  an  oxide 
of  Boron. 
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Table  2  Summary  of  Emissivity  Evolution 


Test  no 

Sample 

Power  [kWj 

Time  [min] 

Tpyrom  [C] 

Tcalib  [C] 

Emiss  [-] 

13 

6 

190 

3 

1600 

1484.89 

0.73 

13 

6 

190 

5 

1620 

1502.6 

0.74 

13 

6 

190 

7 

1625 

1507.1 

0.75 

13 

6 

190 

9 

1630 

1511.5 

0.75 

13 

6 

190 

11 

1625 

1507.1 

0.75 

13 

6 

190 

13 

1635 

1515.9 

0.75 

13 

6 

190 

15 

1630 

1511.5 

0.76 

13 

6 

190 

17 

1620 

1502.6 

0.76 

13 

6 

190 

19 

1635 

1515.9 

0.76 

14 

10 

190 

3 

1620 

1502.6 

0.73 

14 

10 

190 

5 

1620 

1502.6 

0.74 

14 

10 

190 

7 

1620 

1502.6 

0.74 

14 

10 

190 

9 

1625 

1507.1 

0.74 

15 

9 

210 

3 

1685 

1560.3 

0.73 

15 

9 

210 

5 

1690 

1564.7 

0.74 

15 

9 

210 

7 

1675 

1551.4 

0.75 

15 

9 

210 

9 

1695 

1569.2 

0.75 

A  similar  calculation  is  done  using  the  one-colour  and  two-colour  values  of  the 
pyrometer.  Assuming  that  the  measured  intensity  is  proportional  to  the  “representative 
emissivity”  and  the  fourth  order  of  the  temperature,  the  fourth  root  of  the  ratio  of  one  to 
two-colour  temperature  values  are  evaluated.  A  typical  result  is  shown  in  Fig.  10  for 
sample  9  (Test  15).  This  plot  is  the  calculated  emissivity  value  vs  front  face  temperature. 
The  arrows  indicate  the  evolution  of  time.  Initial  data  show  an  emissivity  value  of  0.55  at 
-1500K  at  the  beginning  of  the  test.  As  the  time  passes  the  temperature  increases  to 
1850K  and  at  the  very  end  of  the  test,  the  temperature  decreases  back  to  1500K.  But  at 
this  moment,  the  emissivity  is  calculated  to  be  0.59,  which  is  higher  than  the  initial  value. 
So,  according  to  this  plot,  the  emissivity  value  is  higher  at  the  end  of  the  test  for  the  same 
front  face  temperature.  However,  these  results  should  be  interpreted  with  care  as  there 
might  also  be  some  transient  effects  governing  the  behaviour. 


test  15  sample  9 


Tcolihi.ited  [K] 


Figure  10.  Emissivity  vs  Temperature,  from  pyrometrv  results 
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7  REBUILT  RESULTS  AND  ABACUS  PLOTTING 

A  quick  explanation  about  the  catalycity  determination  is  given  here,  more  details 
could  be  found  in  the  references  [5-8]. 

The  VKI  Boundary  Layer  Code  (P.  Barbante,  Ref  [4])  provides  the  computation  of 
the  heat  flux  at  the  stagnation  point  for  defined  conditions  of  pressure  (pe)  and  temperature 
( Te)  at  the  outer  edge,  catalycities  y,  and  temperature  at  the  wall,  Tw. 

The  functional  dependence  can  be  expressed  as  : 

qw=qw(peJe,Tw,S,^-,-^(^p\e,Ve,/i) 

This  system  can  be  solved  when  the  values  for  the  recombination  probabilities,  y,  of 
the  material  in  a  defined  gas  mixture  are  known:  the  assumption  of  fully  catalytic 
material  (copper)  and  the  use  of  an  effective  recombination  probability  for  air  give  y w=7. 

The  Rebuilding  Code  (ref.  [7]),  implemented  in  Fortran  90,  uses  an  iterative 
procedure  to  solve  this  system:  experimental  and  numerical  fully  catalytic  cold  wall  heat 
fluxes  are  matched  by  iterating  on  the  variables  Te  (which  detennine  He)  and  Ve.  He  and 
Ve  are  respectively  the  enthalpy  and  the  velocity  at  the  edge  of  the  boundary  layer.  The 
important  features  and  assumptions  of  this  Rebuilding  Code  are  given  in  ref.  [7,8],  but  the 
most  important  assumption  is  that  the  edge  conditions  are  local  thermodynamic 
equilibrium. 

The  global  methodology  for  the  enthalpy  rebuilding  is  sketched  in  Fig.  11.  Besides 
the  four  parameters  coming  from  the  ICP  Code,  four  other  experimental  (and 
geometrical)  parameters  are  required  as  input  for  the  Rebuilding  Code.  During  the 
execution,  a  first  choice  for  the  temperature  at  the  boundary  layer  edge  (and  other 
options)  is  asked. 


(n+l) 


He,Vs 


Figure  11.  Enthalpy  rebuilding 
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Note  that  this  Rebuilding  Code  also  takes  the  Barker  effect  into  account,  thus  no 
correction  is  needed  on  the  input  parameter  Aps.  The  Barker  effect  is  known  as  the 
viscous  effect  that  affects  the  total  pressure  measurements  at  low  Reynolds  number  (i.e. 
Re  <100). 


7.1  Rebuilding  Results 

Using  the  experimental  heat  flux,  static  pressure  and  dynamic  pressure  values 
(presented  in  Table  1),  the  boundary  layer  edge  conditions  have  been  calculated  using  the 
methodology  described  in  the  previous  section.  The  results  are  given  in  Table  3,  where 
total  enthalpy,  velocity  and  static  temperature  values  at  the  edge  of  the  boundary  layer  can 
be  seen,  as  well  as  catalycity  values.  Here,  Table  1  is  re-presented  for  the  sake  of  easiness 
in  interpreting  the  results. 


Table  1  Summary  of  Tested  Samples 


Test 

no 

Sample 

Date 

Time 

pstat 

[mbar] 

Power 

[kW] 

H.  Flux 
[kW/m2] 

p  dyn 
[Pa] 

Duration 

Tsteady 

[C] 

5 

20 

31/5/2006 

14:30 

100 

150 

510 

24 

10  mins 

1241 

6 

16 

31/5/2006 

16:00 

100 

160 

715 

25 

10  mins 

1387 

7 

15 

1/6/2006 

10:05 

100 

160 

730 

26 

5  mins 

1356 

8 

14 

1/6/2006 

11:25 

100 

160 

755 

26 

20  mins 

1387 

9 

13 

1/6/2006 

14:00 

100 

170 

905 

28 

10  mins 

1458 

10 

12 

1/6/2006 

15:35 

100 

170 

875 

27 

20  mins 

1436 

11 

11 

2/6/2006 

11:00 

100 

170 

900 

27 

5  mins 

1456 

12 

5 

6/6/2006 

15:45 

100 

160 

710 

24 

16  mins 

1352 

13 

6 

7/6/2006 

11:15 

100 

190 

1095 

31 

20  mins 

1512 

14 

10 

7/6/2006 

14:25 

100 

190 

1050 

31 

10  mins 

1503 

15 

9 

7/6/2006 

16:05 

100 

210 

1225 

38 

10  mins 

1574 

16 

4 

8/6/2006 

12:05 

100 

210 

1165 

39 

15+5  mins 

1547 

Table  3  Summary  of  Rebuilding  Results 


Test  no 

Sample 

Tsteady  [C] 

Qrad  [kW/m2] 

he  [Mj/kg] 

Te  [K] 

Ve  [m/s] 

yH 

5 

20 

1241 

223.522 

14.92 

5437.26 

95.65 

3.227E-04 

6 

16 

1387 

323.023 

20.43 

5872.92 

106.76 

1.351E-03 

7 

15 

1356 

299.563 

20.65 

5887.04 

109.22 

7.577E-04 

8 

14 

1387 

323.023 

21.36 

5932.69 

110.33 

1.108E-03 

9 

13 

1458 

381.929 

25.15 

6159.89 

120.49 

1.483E-03 

10 

12 

1436 

362.882 

24.54 

6124.02 

117.38 

1.256E-03 

11 

11 

1456 

380.168 

25.23 

6164.18 

118.43 

1.465E-03 

12 

5 

1352 

296.632 

20.48 

5876.22 

104.68 

7.891  E-04 

13 

6 

1512 

431.860 

29.85 

6438.52 

134.50 

1.469E-03 

14 

10 

1503 

423.217 

28.60 

6361.95 

132.43 

1.533E-03 

15 

9 

1574 

495.054 

31.90 

6573.08 

152.75 

2.032E-03 

16 

4 

1547 

466.738 

30.09 

6463.69 

151.31 

1.841E-03 
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One  of  the  important  assumptions  in  the  approach  used  for  the  detennination  of 
catalycity  is  that  the  polished  copper  calorimeter  surface  has  y  =  1.  The  true  value  of 
copper  catalytic  efficiency  is  more  likely  to  be  closer  to  0.1 .  To  illustrate  the  sensitivity  of 
the  heat  flux  determination  to  y  =  1  or  0. 1,  the  variation  in  heat  flux  for  varying  catalycity 
is  shown  in  the  following  figure.  The  heat  flux  increases  by  about  2%  for  y  =  1,  compared 
to  y  =  0.1,  and  this  is  already  within  the  uncertainty  of  the  measurements. 


Figure  12.  Heat  flux  as  a  function  of  surface  catalycity. 


8  NUMERICAL  SIMULATIONS 

The  determination  of  catalytic  properties  of  TPS  materials  is  a  major  issue  for  the 
aerospace  vehicles.  These  properties  strongly  affect  the  heat  transfer  to  the  materials  with 
up  to  a  factor  of  two  greater  heat  flux  for  a  fully  catalytic  material  compared  to  a  non- 
catalytic  material.  Knowledge  of  TPS  catalytic  properties  is  extremely  important  for 
designing  aerospace  vehicles  that  have  very  stringent  mass  budget,  for  reusable  launch 
vehicles  the  problem  is  even  more  critical.  Two  main  difficulties  arise  to  carry  out  these 
investigations  and  tests.  Firstly,  the  complete  duplication  of  the  aerothennodynamic 
conditions  for  hypersonic  flight  cannot  be  achieved  in  existing  ground  test  facilities.  Thus, 
only  the  local  heating  environment  on  critical  parts  of  the  spacecraft  can  be  replicated  in 
discharge -plasma  facilities  for  sufficient  test  times  to  evaluate  TPS  material  perfonnance 
in  real  flight  conditions.  Secondly,  the  partial  catalycity  effects  involve  wall  chemistry  in 
non-equilibrium  condition.  This  fact  makes  detailed  measurements  not  easy  to  realize  and 
in  most  of  the  cases  the  investigations  lead  to  a  global  information  for  which  a  data 
processing  is  need  to  determine  the  chemical  contribution  to  the  total  heat  flux  at  the  wall. 
For  this  testing,  a  methodology  centred  on  Local  Heat  Transfer  Simulation  (LHTS),  has 
been  developed.  It  allows  evaluating  the  catalytic  properties  of  TPS  materials  based  on  a 
combination  of  Plasmatron  tests  and  CFD  modelling.  A  fully  catalytic  stagnation  point 
heat  flux  probe  and  a  pitot  probe  are  used  for  the  characterization  of  the  plasma  flow  at 
the  location  of  the  test  sample,  for  given  Plasmatron  conditions.  The  numerical 
computations  of  the  same  test  configuration  provided  the  aerodynamics  parameters  at  the 
edge  of  the  boundary  layer  of  the  sample.  These  two  sets  of  result  are  then  used  as  inputs 
data  for  a  reacting  boundary  layer  code  which  provides  the  heat  flux  abacus  (heat  flux  at 
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the  stagnation  point  versus  wall  temperature)  with  wall  catalycity  coefficient  as  a 
parameter.  The  effective  catalycity  coefficient  of  the  TPS  material  can  be  determined  by 
interpolation,  in  this  reference  abacus,  plotting  the  stagnation  point  heat  flux  and  the  wall 
temperature  obtained  during  the  sample  test  measurements. 

This  report  documents  the  properties  of  the  flow  field  for  two  representative  test 

cases. 


8.1  ICP  Simulation 

For  numerical  investigations  of  the  flow  field,  two  experimental  conditions  have  been 
chosen.  Firstly  those  of  test  no.  9  and  secondly  those  of  test  no  15.  They  are  summarized 
as  follows: 

■  Mass  flow  rate  =  16g/s, 

■  Pstatic  =  lOOmbar, 

■  Pem,9  =  1 70kW  and  Pem,i5  =  2 1  OkW. 

The  ambient  conditions  used  in  the  simulations  of  the  flow  inside  the  PLASMATRON 
coincide  with  the  conditions  actually  used  in  the  experiments  with  the  exception  of  the 
imposed  electromagnetic  power.  The  generator  efficiency  is  generally  assumed  to  be 
-55%. 


The  simulations  were  carried  out  assuming  a  local  thermodynamic  equilibrium  air  model. 
The  transport  properties  are  computed  by  the  help  of  kinetic  theory  and  the 
thermodynamic  properties  by  the  use  of  statistical  thermodynamics.  The  computations  of 
both  sets  of  properties  are  done  by  an  in-house  developed  library  called  PEGASE.  The  air 
model  consist  of  a  mixture  of  eleven  species  as  follows:  CE,  N2,  NO,  O,  N,  N+,  0+,  NO+, 
Ni+,  O2"1",  e  . 

The  Non  Dimensional  Parameters  obtained  are  shown  in  Table  4.  They  serve  as  the 
numerical  input  for  the  boundary  layer  stagnation  line  computation.  This  set  is  used  for  all 
boundary  layer  computations  since  it  has  been  shown,  that  the  impact  of  them  onto  the 
solution  is  little. 


Table  4  Obtained  non-dimensional-parameters  (NDPs) 
A _ j i_u _ Uly _ Yr 

0.4596  0.2041  0.4489  0.2483 


The  results  of  the  simulations  are  shown  in  Figures  1,2  and  3  in  the  Appendix  A  of  the  file 
UHTC_CFD_Report_VKI.pdf  for  test  no.  15.  The  corresponding  tecplot  files  for  both 
cases  are  delivered  in  the  CD  attached  to  this  report. 

8.2  Boundary  Layer  Properties 

The  distribution  of  the  species  along  the  stagnation  line  is  shown  in  Figures  13  and 
14.  Furthermore  the  temperature  is  drawn.  As  expected  a  significant  amount  of  atomic 
nitrogen  and  especially  oxygen  is  present  at  the  wall. 
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Figure  13  Temperature  and  mass  fraction  in  the  boundary  layer;  Test  no  9. 
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Figure  14  Temperature  and  mass  fraction  in  the  boundary  layer;  Test  no  15. 
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To  find  the  number  densities  that  correspond  to  the  species  mass  fractions,  the 
following  relation  is  used 


n.  = 


PNac, 

MW, 


where  n,  is  the  number  density,  NA  is  Avogadro’s  number,  c,  is  the  mass  fraction,  and  MW 
is  molecular  weight  of  species  i.  If  it  is  not  explicitly  given  in  the  output  files  of  the 
simulations,  the  density  can  be  found  by  using  the  chamber  static  pressure,  the  computed 
temperature,  and  the  gas  constant  for  the  mixture  Rm=ZcR[. 


8.3  Files  Delivered  with  the  Report 

Table  5  summarizes  the  list  of  files  delivered  with  this  report. 


Table  5  Overview  of  delivered  files 


name  of  the  file 

format 

comment  on  the  content 

tecplot 

flow  properties  of  PLASMA!  R(  )N  domain  test  no.  9 

tecplot 

flow  properties  of  PLASMATRON  domain  test  no.  15 

tecplot 

flow  properties  stagnation  line  te: 

st  no.  9 

tecplot 

flow  properties  stagnation  line  tes 

t  no.  15 

Excel 

overview  experimental  and  computational  r 

esults  for  all  tests 

pdf 

this  report 

9  CONCLUSION 

The  tests  conducted  in  the  VKI  Plasmatron  were  done  on  UHTC  samples  for 
pressure  and  surface  temperature  conditions  specified  by  the  customers.  Wall  temperature, 
mass  loss,  heat  flux  and  pressure  measurements  are  perfonned  during  these  tests  and 
during  the  facility  calibration.  Enthalpy  rebuilding  is  followed  by  catalycity 
determination. 

LIF  measurements  were  made  for  N  as  the  target  species,  since  226nm  mirrors  were 
delayed  (not  delivered  on  time).  Test  conditions  will  be  re-done  for  oxygen  LIF  and  O 
species  population  and  temperature  at  the  edge  of  the  boundary  layer  will  be  measured  to 
verify  CFD  results. 
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